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Abstract: Motivated by an open question raised in recent times regarding the phase 
transition during the collapse of a neutron star to form a black hole and related stability 
issues, we have constructed charged neutron stars in AdS^ and show that these stars become 
unstable at a particular value of their radius, regarded as the Chandrasekhar radius. We 
reproduced the calculations recently done in [20] in our AdS^ charged star. The analysis 
shows that the non-Fermi liquid behavior found there in AdS^ is still true in this higher 
dimensional case with the presence of Kosevich-Lifshitz oscillations. 
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1. Introduction 

In recent times there has been a flurry of activity in the interplay between condensed 
matter physics and AdS/CFT trying to bridge the gulf that exists in between them giving 
rise to several interesting original ideas in this frontier (For an introduction to this field 
see for exmple [1, 2, 3, 4]). Recently in a series of papers [5, 6] the authors addressed 
the issue of a holographic description of an astrophysical phenomena i.e. the collapse of 
a neutron star toward the formation of a black hole beyond the so-called Oppenheimer- 
Volkoff (OV) limit [7]. Near the very end of their article [5], they have proposed an open 
question regarding the CFT meaning of the OV limit under consideration. The Neutron 
Star to Black Hole formation gives a hunch of a second-order quantum phase transition 
from a condensed matter physicist's view point turning a high density baryonic state into 
a thermal QGP state and thus the issue needs proper investigation as suggested in [5, 6]. 
The main point of investigation is rather open ended and it would be really interesting to 
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settle the issue whether a rehable reahzation of a degenerate fermionic star can indeed be 
achieved conclusively in the AdS/CFT duality framework. We have tried to formulate the 
very problem on our own way rendering hope to give a partial, if not complete, explanation 
of the underlying phenomena which of course needs further study. 

In this note we study charged neutron stars in AdS^ space and focus on their collapse 
toward the formation of a black hole. This process is particularly interesting and worth 
studying from the holographic point of view (see [5, 6] and references therein). This is 
because in terms of the boundary theory this process might as well symbolize a phase 
transition in AdS space. 

First we present a construction of a charged star in AdS space by adopting numerical 
methods. Then we find that such a star ceases to exist at a particular value of the mass, 
charge and radius. This fact implies that this star becomes unstable at this value of the mass 
and charge. This instability will trigger the formation of a new phase and (presumably) 
the star will collapse to form a black hole. Now the behavior of a test fermion in the black 
hole background has recently been analyzed in [8] . It is concluded from their analysis that 
in this black hole phase the boundary theory has excitations about the Fermi-surface which 
is unlike the Fermi liquid. Further research on this point [11] has shown that the Fermi 
vector depends exponentially on the scaling dimension if one uses the duality to extremal 
RN black holes while computing the non-Fermi liquid Green function. This proved the fact 
of violation of Luttinger theorem incase of non-Fermi liquids which is even true in case 
of the extremal AdS^ RN Black Hole. All such computations reported in [8, 10, 12, 14] 
revealed the presence of a Fermi surface in the bulk of a non- Fermi liquid type. For a much 
better and clear understanding of their results and the underlying phenomena that occurs, 
we thought that it might be worthwhile to explore the boundary description of a simpler 
system, that of an ideal Fermi gas in AdS without the presence of a Black Hole. ^ 

We perform a similar kind of analysis by putting in a test fermion in the bulk with 
the charged neutron star serving as the background. From the dynamics of such a bulk 
fermionic field one can compute the two point function of the boundary composite operator 
that corresponds to the bulk fermion. Although, we have not written down the boundary 
correlation function explicitly in the star background, but we gave a possible direction 
how to do that in practice. ^ We hope that, from the pole structure of these two point 
function, it will be evident that in the hydrodynamic limit one gets either a Fermi liquid 
or a non-Fermi liquid type of behavior, nevertheless the explicit computation needs to be 
performed. Via the AdS/CFT dictionary for fermions [9] the conformal dimension of the 
dual operator in the boundary CFT can be controlled by the bulk fermionic mass. In [10], 
it has been already found that by tuning the mass one can match the conformal dimension 
of the boundary operator to that of the free fermion, and thus one obtains the spectral 
function exhibiting a peak which is consistent with a Fermi liquid behavior. It is quite clear 

^It is worth pointing out here that in [19] the authors evaluated the fermionic correlators on a background 
with fermions but without a Black Hole horizon. 

■^The computation is difficult in the sense that the boundary theory under consideration is strongly 
interacting, and the neutron star to black hole collapse process is time-dependent and hence the radial 
holographic dimension can't be easily constructed on the CFT side. 
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from their analysis that one can deviate away from the Fermi hquid behavior by tuning 
the mass away from the free value as the mass was interpreted as the proxy for coupling 
in [10]. 

However, we recalled the fact that if one computes the fermionic correlation functions 
in charged AdS^ black hole backgrounds using probe fermion and analyze the spectral peak 
as done in [11], one finds that the Fermi momentum depends exponentially on the probe 
fermion mass interpreted as a proxy for the coupling and thereby violates Luttinger theorem 
indicating a non- Fermi liquid behavior as was found in [8]. We then show, following [20], 
the presence of Kosevich-Lifshitz (KL) kind of oscillation in our AdS^ charged neutron 
star setting. One must emphasize the fact that , we are working with "global AdS" i.e. 
the spatial sections are rather than planes which implies that the field theory under 
consideration is on time x instead of being on Minkowski space. In [18, 19, 20] the 
authors studied electron stars working with a field theory defined on the Minkowski space 
whereas the authors of [5, 6] adopted our point of view. 

2. Construction of charged neutron stcirs 

In this section we present the construction of the charged neutron star in AdS^. We will 
work with units in which the AdS radius is set to unity. We consider the following metric 
ansatz 

ds'^ = -A{r)dt^ + B{r)dr'^ + r^dul (2.1) 

Besides the metric we have a dynamical gauge field (H) and a fermionic field (ip) in this 
background. We will work in the gauge where the radial component of the gauge field is 
zero (i.e. Hr = 0). With this choice of gauge we write the gauge field as 

H{r) = h{r)dt. (2.2) 

Further we consider fermionic species with mass m, charged under this gauge field with 
charge q. However, we shall treat the fermions in a hydrodynamic approximation. For this 
we will take the limit ^ that the number of particles per AdS radius is infinite with ep/m 
fixed, ep being the Fermi energy. The ideal degenerate Fermi gas formed by these particles 
can then be described by the following hydrodynamic stress tensor for ideal fluids 

T^u = {p + p)utiUi, + pg^u, (2.3) 

where the pressure p and the density p are related to each other by the fiat space equation 
of state for the fermions i.e. 

p + p = fin (2.4) 

'^before taking this limit we must also ensure that the central charge of the boundary CFT is taken to 
infinity (something similar to the large N limit for A/" = 4 Yang-Mills) with all other quantities keeping 
fixed. This is required to ensure that all the multi-trace operators are suppressed. One way to think of 
this limit in the bulk is to consider N species of fermions with the same mass and then take N to infinity 
keeping everything else fixed. 
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The use of the flat space approximation is again justified as there are large number of 
fermions within a given AdS radius and therefore the fermions do not see the curvature of 
the space time. In order to be consistent with this approximation, the mass and charge of 
the fermions must be greater than the AdS radius (i.e. unity in our case). This equation 
of state is implicitly given by 

_ (/X - mf {8m^ + 3fi{3m + /i)) _ 2m^ - Sm^^ + 
P~ 1920^^2 ^P- ' (2-5) 

where /i is the chemical potential for the fermions, which can be expressed in terms of the 
Fermi momentum (kp) by the relation 



// = Y A;|, + m^. 

In this flat space approximation for the fermions the particle number density (n) is given 
in terms of the volume of the Fermi surface and we have 

(2.6 

1287r2 ^ ' 

Then it immediately follows that the charge density (n) is given by 

Note that in our problem these relations are locally valid and generally fi,n,p and p are 
functions of r-coordinate ^ which we ultimately solve for. Here we would like to draw 
attention of the reader to a subtle point, as defined in (2.3) is a static velocity field: 
u^dx^ = A{r)dt and the radial profiles of p and p are determined by imposing the stress 
energy conservation which leads to the following condition: 

dp 1 dA , , , , 

dr A dr 

This equation is very easy to solve. By making use of (2.4) one easily verifies the fact that 
(2.8) is satisfied while the chemical potential obeys 

M(0 = ^ (2.9) 

where at this stage ei;' is an arbitrary constant. Thus the radial dependence of the chemical 
potential pL is simply due to the gravitational redshift. 

2.1 The Equations to be solved 

In the above set up we now write down the dynamical equations which we must solve in 
order to obtain the Neutron star numerically. 



*the symmetries of AdS prevent these quantities from becoming functions of the other coordinates 
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Firstly we have the Einstein equations. The two non-trivial Einstein equations are 
obtained from the rr-component and the tt-component and they are respectively given by 



3r {2A'{r) + rh'{rf) - AA{r) {B{r) {r^p{r) + Gr^ + 3) - S) = 0, 
A{r) {6rB\r) - 4B{rf (rV(r-) - Gr^ - 3) - 12B{r)) - 3r'^B{r)h'{rf = 0, ^'^'^^^ 

where as mentioned before we consider the pressure and the density as functions of the 
radial coordinate. 

Then we consider the Maxwell equations. In this case the non-trivial equation stems 
out from the t-component (which is a mere generalization of Coulombs law) . This equation 
is given by 

A'(r) 6\ , B'(r)h'(r) , . 

^ ' ' h'{r) H ^ ; ^ ' - 2B{r)n{r) - 2h"{r) = 0. (2.11) 



A{r) r J B{r) 

Finally we have to consider the equation of motion for the fermions. However, since the 
fermions are treated in a hydrodynamic approximation this equation is the conservation of 
the stress tensor^. In this case the radial component yields the non-trivial equation and is 
given by 

A{r) {2rB{rfA'{r){p{r) + p(r)) + 'irB' {r)h' {rf + 

-6B{r)h'{r) {rh"{r) + 3/i'(r))) + 3rB{r)A' {r)h' {rf + irA{rfB{rfp{r) = 0. ^'^'^'^^ 

In these equations both p{r) as well as p{r) are present explicitly. However, we will 
eliminate both p{r) and p{r) in terms of /x(r) with the help of the equation of state (2.5). 
Now these equations are extremely non-linear and we have to resort to numerical means 
in order to solve them. 

2.2 Solving the equations 

One obvious solution to these equations is the charged black hole in AdS^ with the chemical 
potential being constant throughout the space. In terms of the above mentioned functions 
this solution may be written as [13], 



Mr) 



B(r)={l+r^l-^ + 9l 



(2.13) 



n{r) = m, 

h{r) = pB-^- 

where m and hb are constants. The parameter m is related to the ADM mass of the Hole 
as M = -^g'^^ m [13] where in natural units IGttG = 1 and oo^ is the volume of the 3-sphere. 
Also, /i_B is the electrostatic potential difference between the horizon and infinity. Since the 
constant value of the chemical potential outside the star is m, the mass of the fermionic 



^which is the relativistic version of the Navier-Stokes equation. 
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species, therefore we use the same notation here also. Similarly we use //_b to denote the 
constant part of the gauge field as it ultimately turns out to be the boundary chemical 
potential. 

Now it is expected that outside the neutron star our solutions should reduce to the 
black hole solution. Therefore, we shall obtain a solution inside the neutron star and then 
patch up our solution with this black hole solution outside. Thus we start with a boundary 
condition at the origin (which may be thought of as the centre of the star) and make a 
choice of the time coordinate inside the star such that at the boundary the first derivative 
of the field strengths match. In order to determine the consistent boundary conditions at 
the origin we solve the equations about r = and find that. 



180 V 87r2 327r2 
90'- ( 16^^-12^ + 30-^°) (2.14) 



r 



2 



'"(^) = + ooi«'n 2 A (32AomVo - 120^omVo + 160^o"iVo " 72Ao/i^ 
- 460807r2AoAio + ISSm'^g^ _ 27Qm'^q^^l + ISSg^^^) + 0{r'^), 

solve the equations (2. 10), (2. 11) and (2.12) upto 0{r'^). Here Aq and /io are the values 
of A{r) and /i(r) at the origin. These are the parameters of our problem and we have to 
choose values for these parameters which serve as initial values of our differential equations. 
Again m and q are the mass and charge of a single species of fermion. The most striking 
thing to note about this solution is that the value of B{r) at the origin is fixed to be 
unity. We do not have the freedom to choose this value on independent grounds. Also the 
parameter q in (2.14) yield the charge of the black hole as Q = ^^g^ 1 [13]. 

2.2.1 Numerical Solution of the equations 

We proceed to solve these equations numerically in the following way. At first we fix the 
value of the chemical potential at the origin to be ^q. Then we fix the value of Aq to be 
unity and then fix a scale for the time coordinate in the patch inside the star so as to meet 
the boundary conditions A{r) = 1/B{r) at the radius of the star (R). The radius of the star 
is obtained from the value of r where the density p{r) goes to zero or the chemical potential 
/i(r) goes to m as can be easily checked from the equations (2.5) and (2.13). Further from 
the value of A{R) (which is the same as 1/B{R)), together with the matching condition 
for the Electric field at r = i?, we determine the mass and charge of the black hole solution 
with which we patch up outside the neutron star. Note that in this procedure, matching 
of the first derivative of A(r) (which is the gravitational field strength in a rough sense) at 
r = i? is automatic. As a part of our choice of units we take the AdS radius to be unity. 
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Figure 1: Solution of various functions inside the star. The values of the parameters are m = 3, 
q = I and /Xq = 6. We have shown the plots upto r=1.48 which is the critical radius of the star. 

In fig:l, we present the solution when the parameters m = 3, q = 1 and fio = 6. This 
corresponds to the core density of 2.97. The Mass of the star is 0.073 and the charge of the 
star is 0.0083 ^.The radius of the star is 1.48. Note that the density function p{r) goes to 
zero at the edge of the star. Also the values of A{r) and 1/B(r) match at the end of the 
star. 

2.2.2 A class of Neutron stars 

Now we can vary the value of the chemical potential at the origin to obtain a large class of 
neutron stars. This amounts to obtaining neutron stars for various values of core density. 
Further under a variation of the core density (which varies form to oo) we can make a plot 
of the mass and radius of the star (see fig:2). In this plot we find there exists a maximum 
value of mass for the star sometimes termed as critical mass of the star. ^ Further the 
solution ceases to exist at a particular value of the mass and radius. This is the signature 
of a critical behavior and occurs when the density at the origin approaches infinity. At this 
critical point the neutron star is expected to start collapsing into a black hole. 

In a similar way, we can obtain a plot for the mass vs charge of the star (see fig:3). In 
this plot we see that for a given value of charge two solutions exist corresponding to two 
different masses. Presumably the one with more mass is the stable one, the other being 
unstable. 

®These values of the charge and mass should not be directly compared with the mass and charge of a 
single fermion since in this case we should also consider the Newton's constant which we have not included 
in the present analysis. 

^We can render a holographic interpretation to this critical mass, too. As found in [6] the limiting mass 
in the boundary CFT theory translates to a limiting conformal dimension of the composite operator made 
out of the fermionic primary fields which in the large N limit construct our model degenerate star. 
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Figure 2: Plot for Mass vs the radius of the neutron star as the core density is varied. Here also 
we have taken m — 3 and q = 1. 



We again remind the readers of the intriguing fact that working in "global AdS" we 
have as spatial section thus providing a scale in the field theory, which is the radius say 
R of S"^. There is yet another scale in the problem, the chemical potential ji. The related 
phase transition from the star to a black hole occurs due to the competition between the 
two scales. In particular at the critical point of phase transition one would get Rga ~ ^. 
So the very existence of Rgs is crucial for our purpose. 

In recent years, all the various works [8, 11, 18, 20] on Fermi surfaces consider Minkowski 
space while considering the boundary field theory. This was important in order to have a 
well-defined momentum k. On a sphere k is not a good quantum number, in particular 
if Rg3 ~ ^ one might worry about whether the notion of a Fermi surface makes sense in 
the momentum space. The answer to this is, yes, indeed the notion can be justified. As 
we mentioned in Sec. 2, that in our choice of units the charge and mass of the fermion are 
large, so that the boundary field theory fermions are locally in a flat space making sense 
of a Fermi momentum in approximate terms. 



3. Holographic dual of the charged star 

We now study the holographic description of this neutron star solution. By studying the 
dynamics of a probe fermion in this background one could get a hint at the boundary 
description of this charged neutron star. 

3.1 Dirac Equation in Spherically symmetric space time 

We first present the Dirac equations in spherically symmetric space times which we subse- 
quently solve in our star background [15]. ^ 

We have to incorporate the interaction of the probe fermion with the fermions forming 
the star. This we accomplish by considering a new gauge (different from the existing one) 

*One should pay attention to the fact that while the fermions are free in tlie bulk, they are not so on 
the boundary (i.e. they don't obey Dirac equation on the boundary) as the boundary theory is strongly 
coupled. 
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Figure 3: Plot for charge vs mass of the neutron star as the core density is varied. Here also we 
have taken m = 3 and q = \ 

with only a non-zero time component being the chemical potential, due to the fermions 
forming the star. We justify this procedure in Appendix A. 
We consider the metric to be of the form (as in (2.1)) 



ds^ = -A{r)dt^ + B(r)dr2 + [dQ'^ + sin^ Qd(^^ + sin^ (/)sin2 Qd)?) 
The natural choice of vielbein basis vectors are 



A\ e^ = V-B; = r; e? = rsin^; e;,^ = r sin sin < 



(3.1) 



(3.2) 



where {t, f, 0, A} are the tangent space coordinates. Then the non zero components of 
the spin connection are 



t.f 


1 A'{r) _ 
2^A B' 


UJg' 


-_ 1 




sin 6 sin (p 




sin 6 


Vb ' 


Vb 


COS 


smcp; LjJ^ = 


cos 4 





(3.3) 



Using the above equations the Dirac equation 

{P - m)^} = 0, 

reduces to 

1 1 [A'{r) 3 



cos ( 



r^dr + T'^Vt - ( ^ + - ) r + 

B VA 2Jb\2A r I r 



m h/; = 0. (3.4) 



where ^ is the Dirac operator on the three sphere. The eigenvalues of this operator 
are [16] 

(3.5) 



K = ±i\^-+nj, n G {0,1,2,...}. 

We shall denote the eigenfunction corresponding to these eigenvalues by f^- Also the Vf 
operator denote the gauge covariant derivative and contains both the gauge field and its 
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time derivative. All the gamma matrices in (3.4) have a tangent space index and hence are 
flat space gamma matrices. Here we use the following basis for the gamma matrices. 

r = I®(-a3); r*=I®(zai); = a' a2. (3.6) 

where cr* are the ordinary Pauli matrices. Further we consider the following separation of 
variables 

V' = /. (^^^^e—*. (3.7) 

where (j)± are the eigenvectors of with eigenvalues ±1. Using this separation of variables 
the Dirac equation (3.4) reduces to the following set of coupled first order equations 



^ ifA'ir) 3\ r^fuir) k 



(3.8) 



where k = in, and u{r) = uj + n{r) — 3/2qh{r) ^. 
3.2 The flow equation 

In this subsection we will derive and analyze the so-called flow equation for the Dirac 
fermions in AdS^. Using the set of equations 3.8, one can define the wave function at small 
r to be of the form: (Note that m is the fermionic mass, and it enters in the expansion) 



Q r 



mL , \ _i_ X ( JmL-l) 



7(r~™^ + ...) +5 r(™^-^) + 



(3.9) 

Here L is the AdS radius that we have taken to be unity L = 1. The coefficients a, f3, 7 
and S are what the authors of [11] call A, B, C, D in their article and they are related with 
one another. The retarded Green function could be written as [8, 11] 

Gr = e-2'^lim " 1 (3.10) 




where ^_ is defined as suitably defined ratio between 

= -i^ (3.11) 

From the system of equations, (3.8) the flow equation can be derived, dividing the second 
equation by and inserting into it the first one, finally yielding: 

- 2«_»^/B .e^m-^)^-^{^^'-)^<^ (3.12) 



The numerical solution of this equation has some singularity problems in r, due to the 
choice of boundary conditions, thus preventing its use in the calculation of the correlator 
in the neutron star background. This issue needs further careful investigation. 

®Note that the unconventional factor of 3/2 is present because here we are using a different normalization 
of charge. 
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Figure 4: Plot of 0+ in the star background constructed in the previous section. The vertical line 
denotes the value of r for the edge of the star. Here the values of the parameters are m = 3, q = 1, 
fiB — 1, k — 3/2, and a; = 3. In the right image we focus on the edge where the first derivative of 
the function is continuous. 

3.2.1 The initial conditions 

Now we have to specify the boundary conditions for the equations (3.8). This is done by 
demanding regularity of the solution near the origin (r = 0). This regularity criterion in 
general depends on the value of k once we fix the mass and charge of the fermion. The 
lowest positive value is A; = 3/2, in which we will focus for our present purpose. In this 
case the regularity at the origin demands, that if (p+ is 1 at the origin then should be 
-1 i.e. 

0+(r = O) = l ; (/._(r = 0) = -l 
We shall use this boundary condition to solve the equations (3.8). 

3.2.2 Numerical solution of the Dirac equation 

The solution of the probe fermion in the neutron star background is obtained numerically 
for m = 3, q = 1 and k = 3/2. It is shown in figs: 4 & 5. Here we use the boundary 
conditions as discussed in §§3.2.1. We at first let the differential equation evolve to obtain 
a solution up to the end of the star and subsequently use the value of (p± at this edge as 
the initial value for the subsequent evolution outside the star. The first derivative match 
of the solution in these two patches is automatic (see figs:4 &: 5) and directly follows from 
the fact that all the functions that appear in the equations (3.8) are continuous at the 
edge of the star. One possible direction, one can take from here, is that, one might try 
to compute the boundary correlation function in the presence of the star and understand 
the behavior of the bulk Fermi surface in terms of the boundary correlators, taking the 
same route outlined in the formalism of [17]. A noteworthy point here, is that the only 
difference, we should care about, is that we need to evaluate the fermionic correlators in 
the neutron star background, whereas there exist large number of literature doing the same 
computation in the background of a charged black hole (See for example [8, 10, 11, 12, 14].) 

3.3 Kosevich-Lifshitz oscillations in AdS^ neutron star 

It has been shown recently in [20] that, the charged star shares some important features 
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Figure 5: Plot of 4>- in the star background constructed in the previous section. Here the values 
of the parameters are also m = 3, g = 1, ^lb = 1, k — 3/2, and w = 3. Again the first derivative of 
<j)- is continuous as can be seen from the plot on the left. 



with the Fermi liquid. In particular, the presence of the Kosevich-Lifshitz oscillations [18] 
in magnetic field has been shown in the case of AdS4 charged fermion star. As said in [20] 
the arguments used to demonstrate this feature could be generalized to higher dimensions, 
showing that even in AdSs Kosevich-Lifshitz oscillations are present, nevertheless the Lut- 
tinger theorem is violated. In order to calculate the magnetic susceptibility x = ~fi|^ 
can use the relation between the free energy and the total charge in the dual CFT as in 
[19, 20]: 

n(xQ= y/^r^a{r)dr (3.13) 
Jo 

where in general the charge density is 

a= r g{E)iy{E)dE (3.14) 

J m 

and the density of states in the full theory is 

g{E) = I3E[E^ - - IB) (3.15) 

The easiest way to generalize the argument as presented in [20] to the higher dimensional 
case is to add an extra dimension, i.e. by considering the cylindrical-like metric 

ds^ = -fdt^ + jdr'^ + r^du}2 + dz^ (3.16) 

and then the integral expressing the charge density is generalized in AdS^ just by adding 
an integration over pz 



-oo J J ~oo J ]^ _j_ g T 



The energy, as in the AdS^ case, is 

Eil,p,p,) = y/p2+p2+£Bioc (3.18) 
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The local physical quantities are related to the dual field theory quantities as [20]: 

J, (3.19) 



Bloc (X \b 



Tloc 



'9tt 



In the following calculation we will drop out this subscript, and restore the correct depen- 
dence on physical quantities back to position at a later point. The Poisson resummation 
over Landau levels gives 

^ = V^/ ^^^/ (3-20) 

Now expanding the Fermi distribution in terms of Matsubara frequencies and transforming 
the integral over i into an integral over energies the integral (3.20) becomes 

^ = .^y dp J dp ^(P^+P^ ^ e__ 

2 ^J-oo J-oo J^p^+pl+m^ E - p{l + i^2T:{n + ^)) 

(3.21) 

This integral could be performed rotating the integration path so that the exponential in 
the numerator becomes real and decreasing, i.e. considering the rotation of the energy 
as E' — ?• Ee^"^/^ if > or ii^ — ?• Ee~™/^ if A; < 0, and remembering that the lower 
integration limit is set to -s/p^ + p1 + m? . In this analytical continuation one has to take 
into account the role of the poles given by Matsubara frequencies. These poles contribute 
only if p^ + pI < ^i'^ — m? and in addition n + l/2>0ifA;>0 and n + l/2<0iffc<0. 
The contribution of the poles is the one giving the oscillatory part of the charge density 
a osc: and thus it' s the only one we will consider next, hence through the Residue theorem 
then the Gosc is 



fc,n>0 



D 



(3.22) 



while the last integration is made in the region D, being the circle of radius \/ p? — vn? in 
the plane {p,Pz)- This expression could be inserted in (3.13), thus obtaining the oscillating 
part of the free energy 



^osc oc Qosc = 2nl3Tp Im V H drr^ V^e*"'=(^'-'"')/^e-^-'K"-|)^A^/s x 

Jo 



k,n>0 



(3.23) 



dpzdp e 



-iTTkip'^+pD/B 



D 



Since fi ^ B, a saddle point approximation could be done. The r dependence of the fast 
oscillating exponential is 

^i.ki,^~m^)/B ^ .^^ expiA;^ (3.24) 
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With this approximation the integral (3.13) becomes 



-me 



('■') 



A;,n>0 



(3.25) 



X / ^ dpzdp e 



where, as usual, r* is the value of the radius where the argument of the exponential has an 
extremum. The last integration could be done passing to circular coordinates in the {p,Pz) 
plane inserting the radial coordinate p = Y^P^~+~pf 



dpzdp e 



~iTTk{p'^+pl)/B 



D 



pdpe 



-iirkp'^/B 



e 2B 



2B 



vr 



(3.26) 



where b = \/ p^ — vrP' has been defined. The oscillatory part of the total charge is then 
obtained introducing this result in (3.25). The summation over the index n could be done as 
usual, giving an hyperbolic sine, moreover the sum over k could be suppressed considering 
only the term /c = 1, since the other terms are negligible. In the end after all dusts get 
settled the oscillating free energy is 



oc Qosc = 27r/3 Im 



m 



.^3/2 



in^i^ ■62^sinf|^ 



Fin .^""^^/^...Im 



Fir*) 



y/A"{r*) sinh 27rm*T/B 

bI 2Trm*T/B 
y/A"ir*) sinh 2TTm*T/B 



n>0 

in/A in^^ -i^^ 

e ' e B e 2B 



g-47r2(n-i)m*T/B 

nAir*) 



sm ■ 



2B 



V2 fnAir*) vr 

cos 

2 \ B 4 



(3.27) 



The function Fir*) contains all the factors coming from the integration over the radial 
coordinates, like \/gir*) and the numerical prefactors. As can be seen, the expression found 
agrees with the standard result for three dimensional systems, apart from an additional 
and unimportant offset of the oscillations. In particular the usual scaling dimension with 
the magnetic field is found, i.e. Cl oc B^^"^ [21] and the same frequency dependence of 
oscillations on the extremal section of the Fermi volume is observed. 



4. Discussions 



In a nutshell the findings of our paper is as follows: We built, using numerical means, a 
charged neutron star in the AdS^ and analyzed in detail the related gravitational collapse 
to form a black hole. We depicted the critical behavior of the degenerate star with few plots 
and then solved the Dirac equation in the spherically symmetric star geometry. We wrote 
down a flow equation governing the dynamics of fermions constructing the charged neutron 
star in AdS^. We encountered few problems while computing the boundary correlation 
function in the presence of the star, we list the difficulties that arise in finding the fermionic 
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correlator in the star background. In a sense, according to [5, 6], it is an open question to 
find a reliable realization of the degenerate star using AdS/CFT duality and find a CFT 
interpretation of the OV limit [7] in the context of collapse phenomena. Here we note that, 
if one computes the retarded Green function 3.10 in the charged AdS^ extremal RN Black 
Hole background then, according to [11], 

kp kF{0)exp{-amL) (4.1) 

i.e. one observes that the Fermi momentum fall off exponentially from the value it takes at 
zero fermion mass with kp{0) ~ 0.8155 and a ~ 0.80 for AdS^ geometry, clearly indicating a 
violation of Luttinger's theorem and hence proving a non- Fermi liquid behavior in the black 
hole phase. This initially prompted us to make a conclusion that one possible holographic 
interpretation one could dub to this collapse process is that when a neutron star collapse 
to form a black hole, seeing from a boundary point of view one finds a second order 
phase transition from an ideal degenerate Fermi gas to a non-Fermi liquid. Recently in 
a paper [20] the authors have found that the Kosevich-Lifshitz oscillations still persist 
in the charged neutron star phase. In the same article the authors also explain how to 
reconcile the violation of Luttinger theorem in the case of the charged fermion star and 
the existence of the KL oscillation and the argument is still true in our AdS^ case : only 
fermions in the spherical shell of radius r* contribute to the oscillations, and thus it' s 
not possible to reconstruct the whole Fermi volume from the analysis of oscillations. From 
the boundary field theory point of view this means that not all the degrees of freedom 
are taken into account through quantum oscillations, and this results in the violation of 
Luttinger theorem. Reproducing the calculations done in [20] in our AdS^ charged star, 
we have found again the appearance of the KL oscillations due to magnetic field consistent 
with the standard results for three dimensional systems (see [21]). After all these, it is still 
not clear how to realize a holographic dual of the AdS^ star, although it seems to us that 
both the star and the black hole are non- Fermi liquid states of matter. 

Finally, despite trying to give a concrete holographic picture of this underlying phe- 
nomena we tried to rethink from a different perspective, all the existing results on this 
subject and did a coherent study of all of them reaching a conclusion that this elusive 
phenomena still require further studies. 
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A. The zero temperature and finite chemical potential two point function 
in flat space 

In our analysis above we have used the crucial fact that the effect of finite chemical potential 
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at zero temperature is captured by introducing a gauge field whose time component is the 
chemical potential (all other components being zero). In this section we shall try to justify 
this statement by considering fermions in flat space with a constant potential. We shall 
do this by considering the two point function of the fermions. At first we shall view, the 
introduction of the chemical potential, as a redefinition of the vacuum state of theory of 
free fermions. Them from there we shall demonstrate that the two point function computed 
in this new vacuum state is the same as that computed for the case where the chemical 
potential is introduced through the time component of the gauge field. The later method 
more elegant and easy to generalize. In fact, we have used this above in more complicated 
settings where the chemical potential has a spatial variation. Therefore, in order to be 
certain, we present an analysis of the situation in this simple setting and verify the 
equivalence of the two methods. 

A.l The free Fermi sea: The operator calculation 

In this subsection we consider free massive fermions in flat space with the fermions being 
filled upto the fermi level with fermi momentum kp. Thus the chemical potential of the 
system is non-zero and is given by 



with m being the mass of the fermion. We consider the system to be at zero temperature. 
We will be interested to calculate the Feynman propagator or the time ordered two point 
correlator or this system. We shall perform this through an operator calculation in which 
the two point function is the expectation value of a product of two field operators in a state. 
This state is the one in which all the single particle states upto the Fermi momentum is 
filled up. 

We consider free fermions in flat space with the system being described by the Dirac 
Lagrangian [22] This theory is a quadratic theory and can be solved exactly. The fields ip 
and ip can be expanded in terms of the creation and annihilation operators as follows 

" / 7^71^3 exp i-ip.x) + bpv'{p) exp (ip.x)^ , 

I (A.l) 



Now we shall incorporate the presence of a finite chemical potential by constructing 
a state in which the fermions are filled upto the energy equal to the chemical potential. 



^°This fact is certainly true and is very well known for finite temperature. In case of finite temperature 
the chemical potential can be introduced but putting a twisted boundary condition for the fermion (instead 
of a mere anti-periodic one) along the compactified time circle. It is very well known that such twist can be 
undone by a gauge field whose time component is the chemical potential. This implies that the chemical 
potential can also introduced through such a gauge field without putting the twist. What we verify here is 
that it this fact continues to hold even at zero temperature. 

^^Here we shall define the chemical potential to be the energy required to add one more particle at the 
fermi momentum 
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characterized by the fermi momentum. Let us denote the state in which the fermions are 
filled upto the fermi level by \kp) such that, 

kp 

IM = n^2£pa;t|0). (A.2) 

p 

Now for the time ordered correlation function we are required to calculate the quantities 

{kF\'ilj{x)'4){y)\kF) and {kF\'4^{y)'4'{x)\kF) ■ These quantities evaluate to 

{kF\i){x)i){y)\kF) = / 77^77Tr(i' + "i)exp(-p.(2;-y)), 
f°° d^p 1 

{kF\^{y)i){x)\kF) = 777^7TTr(i'-"i)exp(-p.(y-x)) 

Jkp {2-KY2Ep (A.3) 

+{p - m) exp {-p.{y - x))) . 
Now the time-ordered two point correlation function is given by 

5-^(3; _ - J i^F\ipi^)'^iy)\f^F), for xo > yo, (close the contour below), 
I — {kF\'ip{y)ip{x)\kF) , for yo > xq, (close the contour above). 

Now the above two point function can be captured in the contour integral 

(27r)4 p2_^2 exp-iff.(x-7/). (A.5) 

where the contour for po is chosen as shown in fig:6. Nevertheless, the same answer can be 
obtained with the usual contour prescription (the dotted red line in fig:6) if we include a 
real shift of fi in po- Therefore it is convenient to define the shifted variable, 

Po = Po - At 

and make this substitution in (A. 4). Then in terms of this shifted variable po the contour 
prescription is the usual one. Now in terms of this shifted variables, the effect of ^ is 
completely captured is we introduce a gauge field whose time component is /U. This is true 
if we consider the eigenvalues of the dt operator to be instead of po- This justifies our 
use chemical potential as the time component of the gauge field. 
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